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ABSTRACT 


A two-layer, nonlinear, equatorial 6-plane model, in 
p-coordinates is used to study the atmospheric response to a large scale 
prescribed heat source varying in time. The heat source is meant to represen 
a convective burst with total duration of approximately 48 hours over the 
Amazon/Bolivia region. The boundary conditions used are meridional velocity 
zero at 60°N and 60°S, oj = 0 at the top and zero geometric velocity at the 
lower boundary. Sensitivity study was done which includes initial state at 
rest. compared with realistic initial flow. The scale of the heat source is 
1500km in latitude and longitude and it is centered at 10°S. Special 
attention is paid to the distribution and intensity of the induced vertical 
motion. 


The model is integrated for two days and the preliminary results 
show agreement with the observed 200mb flow. Of interest is the establishment 
of a trough and descending motion to the northeast of the heat source. 

A conjecture is thus made that the Amazon heat source and its fluctuations 
bear some relationship with the drought problem over Northeast Brazil. 


1. Introduction 


The Northeast Brazil climate has been the subject of study by 
scientists, who have been trying to explain the observed climatic anomalies 
due to the disruption of the lives of over 30,000,000 people. Two main 
lines of research have been followed. The first one considers possible 
local effects as the main cause of the anomaly. Along this line. Gomes 
Filho [4] studied two mechanisms of formation and maintenance of the 
semiarid region: albedo gradients and orographic effects. However this 
conclusion is not well established yet. Along the second line of research, 
the main cause is due to external mechanisms such as those invoked by 
Namias [16], Hastenrath and Heller [8], Kousky [11], Moura and Shukla 
[15] and Buchmann [2]. Namias [16] demonstrated a connection between the 
700mb circulation pattern over the North Atlantic and the rainfall at 
Quixeramobim (state of Ceara). Hastenrath and Heller [8] used a large 
number of observations to show that rainfall in the region of Ceara is 
closely linked to the meridional displacement of the equatorial trough 
zone. Kousky- [-11]' showed that: a) frontal systems penetrate the southern 

part of Northeast Brazil throughout the year; b) frontal incursion 
plays an important role in the December-January maximum in the month 
precipitation in the southern part Northeast Brazil; c) in certain 
periods, the frontal incursions affect rainfall as far north as Ceara. 

Moura and Shukla [15] proposed that a possible mechanism for the occurrence 
of severe droughts over Northeast Brazil is the establishment of a 
thermally direct local circulation which has the ascending branch at about 
10°N and the descending one over Northeast Brazil and the adjoining oceanic 
region. Buchmann [2] studied the occurrence of drought and rain 
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anomalies in the Northeast Brazil, in relation to variations in the 
synoptic pressure systems located in middle latitudes of the Northern 
Hemisphere, through the mechanism of lateral forcing. 

However, there are other mechanisms not yet fully explored. 
Because the amount of latent heat released during the rainy season of the 
Amazon region is so large (Kasahara and Mizzi,[9] ), one can expect a 
significant effect of the local forced circulation on the neighboring 
regions, and assume a direct thermal circulation established over the 
tropical sector of South America, with maximum upward movement over the 
Amazon region and the downward branch of the resulting Walker circulation 
affecting Northeast Brazil and the adjacent Atlantic Ocean (Gill, [3] ). 
This hypothesis agrees with Bjerknes Cl] and Kidson [10] ideas 
concerning the relationship between the upper tropospheric high (also 
known as Bolivian High after Gutmann and Scwerdtfeger , [5] ), and the 
variability of precipitation over Northeast Brazil. Specifically, the 
circulation pattern over South America shows large variations in different 
time scales: (i) from the Southern Hemisphere summer to winter, the 200mb 
circulation changes from a well-defined upper tropospheric anticyclone to 
near equatorial westerlies (Kruels et al.,[13]); (ii) convective bursts 
during summer make a substantial rearrangement of the upper trospospheric 
anticyclone and the accompanying through over Northeast Brazil (Silva Dias 
et al.,[19]; Vi rgi , [21 ] ) , and ( i i i ) the large diurnal variation of 
convection seems to have a strong signal in the divergent component of the 
upper tropospheric circulation (Silva Dias et al.,[20]). Latent heat 
release is a plausible controlling mechanism on the three time scales of 
the tropical tropospheric circulation over South America. Due to the 
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confinement of the precipitation pattern in the central and western part 
of Brazil during the summer months, it is reasonable to investigate the 
effect of latent heat-induced circulations on Northeast Brazil. 

Fig. la shows a typical infrared satellite image with the well- 
developed convection over the Amazon region. Fig. 1b shows the associated 
200mb streamlines and the vertical pressure velocity field (u) at 500mb 
(thin solid and dashed lines for upward and downward motion, respectively), 
as obtained from the National Meteorological Center (NMC) analysis for 
January 5, 1981. Similar upper tropospheric situations as well as the 
transient behavior are shown in Fig. 6 of Silva Dias et al. [19] . 

Upward vertical motion is observed in association with the: upper 
tropospheric circulation and sinking motion occurs over Northeast Brazil 
when the upper circulation is cyclonic. The satellite imagery qualitatively 
agrees fairly well with the NMC vertical motion. In particular, recent model 
studies by Silva Dias et al.[19] indicated that transient convection on 
the typical time and spatial scales over the Amazon region tend to generate 
a Walker- type circulation associated with Kelvin waves to the east of the 
heat source, an upper tropospheric anticyclone centered to the southwest 
of the source and an upper trough over Northeastern Brazil in agreement 
with the observed pattern shown in Fig. 1.. However, the model results were 
based on the linearized shallow water equations about a basic state at 
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rest. The objective of this paper is to further explore the response of 
the atmosphere to transient forcing, in terms of the associated vertical 
motion and possible nonlinear effects with a two-level primitive equation 
model. The time and horizontal scales of convective bursts, centered at 
about 10°S, are assumed to be of approximately 48 hours and of the order of 
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1500 km, respectively. 
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2. Observational Evidence 

The summer upper level circulation over tropical South America is 
characterized by the anticyclonic flow with large transient variations 
(Virji, [21]). In order to analyse the transient behaviour of the upper 
circulation and its association with cloudiness, we obtained the vertical 
motion at 500mb at 00:00 GMT and 12:00 GMT for January, 1981, based on the 
NMC analyses. This choice was based on the particularly intense convective 
activity observed in the period (Silva Dias et a!., [20]). A mass balanced 
kinematic method between lOOOmb and 200mb is applied to the gridded data 
(horizontal resolution of approximately 5° x 5°). Infrared geostationary 
images processed by the Instituto de Pesquisas Espaciais - INPE (Institute 
for Space Research) were collected for January, 1981 at times as close as 
possible to 00:00 GMT and 12:00 GMT. 

The diurnal signal in the vertical motion is clearly evident 
during the period with maximum upward motion at 00:00 GMT over the Amazon 
region and Central Brazil (approximately 20:00 LT). The diurnal variation 
in the vertical motion field ( u> in pressure vertical coordinate) over 
tropical South America has been studied by Silva Dias et al. [20], who 
showed that the signal is significant from a statistical point of view. 
However, as noticed by Virji [21], convective bursts often occur over tropical 
South America, frequently in response to penetrating mid latitude systems 
(Kousky and Virji, [22]). One such sequence of satellite images (Fig. 2) was 
observed from the 22nd to the 25th of January, 1981 with enhanced convective 
activity on the 23rd, as shown in the infrared satellite image at 21:16 GMT 
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(Fig. 2.b). Fig. 3 shows the accompanying upper level wind vectors 
(200mb) and the 500mb vertical motion (solid lines for w £ 0 and dashed 
lines for to < 0). A comparison between the computed vertical motion and 
convective activity (Fig. 3 and 2 respectively) shows a reasonable 
agreement, with maximum upward vertical motion at 00 GMT on the 24th 
(Fig. 3c for the co-field and Fig. 2b r o r the satellite image). 

Several interesting points deserve attention in Fig. 3: (i) 
before the convective burst, the anticyclonic circulation is zonally 
elongated towards the west with center at approximately 18°S and 75°W 
(Fig. 3a); (ii) as the convective activity enhances, a strong cross 
equatorial flow over the northern Amazon region is observed (Fig. 3c), 
accompanied by enhanced upward vertical motion over Central Brazil (10°S, 
55°W) ; (iii) as the convective activity decreases, the easterly component 
is established just off the equator (Fig. 3f); (iv) the center of the 
anticyclonic circulation is displaced eastward as the convective activity 
enhances (Fig. 3a to 3c), presenting a westward elongation after the 
maximum activity (Fig. 3e and 3 f); (v) ridging is observed (25°S, 55°W) 
during the final period (Fig. 3e and 3f) just south of the most active 
region. These general observations were also detected in other periods 
of January 1981 and during a particular period of February 1979 (8th to 
12th) as discussed by Silva Dias et al. [19] . 

3. Governing Equations 

The fully nonlinear primitive equations above an equatorial 
B-plane are discretized in the vertical assuming a two-level model , with 
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the zonal and meridional wind components and potential temperature defined 
at 750mb and 250mb. The geopotential and vertical pressure velocity are 
defined at the bottom pressure level ( lOOOmb) , 500mb and top (Omb). The 
vertical boundary conditions are oj = 0 at the top and zero geometrical 
velocity at the lower boundary (W = 0). Mo topographic effects are included. 
Thus, there are two vertical modes in the linearized version of this model 
(about a basic state at rest): the external (barotropic) and the internal 
mode with reversing sign in the vertical (baroclinic). Schuman's horizontal 
discretization scheme is applied to the governing equations (Haltiner, 

[7] ) and the time differencing scheme is centered in time. The horizontal 

boundary conditions are constant in time with v = 0 at 60°N and 60°S. Zonal 
periodicity at the equatorial circumference of the earth is also assumed. 

The horizontal boundary conditions are not expected to influence the model 
results because the time scale of the integration is sufficiently short to 
avoid boundary contamination. 

The model circulation is forced by a known heat source which is 
Gaussian-shaped in the horizontal direction and with time dependence as 
shown in Fig. 4, with maximum heating at 24 hours. In the particular case 
of the two level model, the heat source is assumed to have the same 
intensity at 750mb and 250mb, corresponding to 6°C/day. The choice of 
the spatial and temporal scales are based on the satellite imagery. The 
heat source intensity is, however, an estimate based on observations in 
other tropical regions (Riehl , [17] ) because there are no aerological 
estimates of the transient characteristics of the Amazon heat source, to 
the authors' knowledge. 
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The results shown in this paper are confined to the model outputs 
at 24 and 48 hours. No initial instabilities were observed in the model 
integration because the transient forcing is slowly increasing in time 
rather than a switch-on forcing. A carefull observation of the vertical 
motion field at frequent intervals (4 hours) shows a very stable and 
continuous evolution of the numerical solution. 

Numerical solutions including Rayleigh damping and Newtonian 
cooling were obtained using a characteristic time scale of 10 days. Although 
- the flow intensity was somewhat reduced compared to the frictionless version, 
no major discrepancies were observed up to 2 days of integration time. 
Constant radiative cooling of 1.5°C/day was also tested in the model. Since 
no available potential energy is generated in this case and the model 
contains no moisture physics or radiation interaction, no horizontal motion 
is generated by such forcing. In fact, the prescribed heat source can be 
interpreted as the combination of both latent heat release and radiative 
vertical flux divergence is certainly an important process (Gray and 
Jacobson, [6]) due to the typical difference between the vertical profiles of 
clear air/cloudy stay radiative cooling rates. In order to keep the thermal 
balance, clear regions are associated with sinking motion which approximately 
balances the radiative cooling. Horizontal motion is also generated by the 
difference between the radiative cooling rates in cloudy and clear regions, 
but this effect is included in the prescribed heat source. Thus, the 
results shown below correspond to the frictionless adiabatic version. 

4. Model Results 

Silva Dias et al. [19] computed how the forced energy is 
partitioned between Kelvin, mixed Ross by- gravity, Rossby and/gravity modes 

/ 



- 8 - 

for the linearized primitive equations about a basic state at rest 
assuming different temporal and spatial scales for the heat source. The 
characteristics of the response are quite different from the stationary 
forcing (Gill, [3] ; Moura and Shukla, [13]). The main findings are: (i) 
Rossby wave energy increases as the spatial scale and the latitude of the 
forcing increases; (ii) as the forcing time scale decreases, less energy 
is excited in the high frequency modes; (iii) as the forcing spatial scale 
increases, the energy in gravity modes decreases sharply; (iv) for 
fixed time scale of the forcing, the energy in high frequency modes 
initially increases with the horizontal scale of the forcing and then 
decreases. 


The results obtained with the two-level nonlinear primitive 
equation model are essentially similar to the linear situation shown in 
Silva Dias et al . [19] , primarily concerning the rotational contribution 

and the Kelyin component to the east of the source. Figs. 5 and 6 show the 
250mb wind field at 24 and 48 hours, respectively , for a heat source 
centered at 10°S, with symmetrical e-folding width of the order of 1500km. 

The corresponding vertical motion field at 500mb is shown in Figs 7 and 8. 

The development of the upper tropospheric anticyclone and the trough over 
Northeast Brazil is well depicted by the model (compare with Fig. 4 of 
Silva Dias et al.H9] ). The heating over the Amazon region initially 
induces subsidence over Northeast Brazil, although the forcing function 
is quite wide. It is interesting to note the preferred location of the 

i 

downward motion. As time increases, a slight ascending motion can be j 

detected over Northeast Brazil. However, with smaller scale forcing, the j 

/ 

subsidence over this region is emphasized as shown in the vertical motion 
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field at 24 hours for a heat source located at 15°S and with a 
meridionally elongated source (zonal and meridional e-folding width of 
700 kn and 1500 km, respectively), as shown in Fig. 9. 

Possible nonlinear effects were also studied by prescribing an 
initial condition with a more realistic zonal flow with subtropical jets 
located at about 35° to the north and south of the equator, attaining 
35 ms -1 . In this case the wind field at 48 hours at 250mb is shown in Fig. 
10 and the deviation from the initial flow is shown in Fig. 11. The basic 
characteristics of the forcing are the same of Figure 6. The main 
difference from the no initial flow case is due to the stronger meridional 
component to the south of the upper anticyclone. This effect implies a 
larger distortion of the zonal flow (ridging) and therefore we might 
expect important downstream amplification of the wave train if longer 
time scale forcing is applied or if a succession of convective bursts occurs 
(Lau and Lim, [M] ). The upper westerlies to the east of the heat source 
are also evident in Figs, 5, 6, 9 and 10. This upper level configuration 
seems to be related to free Kelvin modes according to Silva Dias et al., 
[19] an d recently found in satellite data by Salby et al [18] . 

The model results seem to corroborate observed characteristics 
of the tropospheric circulation over tropical South America during 
actively convective episodes, as presented in Section 4. In particular, 
there is correspondence with the intense cross equatorial flow at the 
initial stages of the forcing (Fig. 5), the slow turning of the wind to 
easterlies south of the equator (Fig. 6) and the westward displacement of 
the upper anticyclone center at the latter stages (Fig. 6). The case with 
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initial zonal flow indicates that nonlinear interactions might be relevant 
in emphasizing the observed ridging at the final stages of the convective 
burst, as observed in Figs. 3e e 3f. 

5. Conclusion 

The nonlinear primitive equation model simulation of the transient 
heat source reproduced most of the results obtained with the linear shallow 
water equations of Silva Dias et-al. [19] . However, the results also 
indicate that the induced subsidence occurs preferentially over Northeast 
Brazil primarily for meridionally elongated heat sources. These results 
indicate that a regional Walker type circulation with ascending motion 
over Northeast Brazil may take place. Thus, excessive convection over 
central and western South America may induce an unfavourable dynamical 
situation for convection over the Northeast region. The model is able to 
simulate the basic features of the observed characteristics of the 
transient convection episodes (Section 2). The nonlinear effects up to 2 
days seem to be confined to higher latitudes where the upper anticyclonic 
circulation interacts with the upper westerlies, intensifying the 
ridging mechanism. The intensification of the equatorial westerlies to the 
east of the forcing (Fig. 6)doesnot seem to be evident in the observations 
(Fig. 3) although westerlies are observed. This can be due to the 
rotational constraints imposed on the analyses and assimilation scheme 
used by NMC. It is interesting to note that transient forcing is more 
realistic than stationary forcing over tropical South America. According 
to the linear model partition of energy, discussed in the previous 
section, more energy in fast modes is expected for highly transient 
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forcing and therefore more intense vertical motion teleconnections are 
expected through gravity wave activity. 

Similar model simulations have also been performed with zonally 
elongated heat sources located over the eastern equatorial Pacific Ocean. 
These model results are expected to represent the El Nino type of anomaly. 
The model results also seem to indicate a preferred subsidence region 
over most of equatorial South America on the time scale of two days 
(Fig. 12). These results are in agreement with the observed negative 
correlation between precipitation over the eastern Pacific Ocean and over 
tropical South America (Kousky et al., [IE]). 
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Figure Legends 


Fig. 1 
Fig. 1 

Fig. 2 


Fig. 3 


Fig. 4 
Fig. 5 

Fig. 6 
Fig. 7 
Fig. 8 
Fig. 9 



- Infrared satellite image for 01/05/81 at 21:16 GMT. (SMS-2). 

- 200mb streamlines (solid lines with arrows) and 500mb vertical 

_4 “1 

motion in 10 mbs on 01/06/81 at 00: GMT. 

- Infrared satellite image of (a) 01/23/81 at 12:16 GMT (SMS-2), 

(b) 01/23/81 at 21:16 GMT, (c) 01/24/81 at 12:16 GMT, (d) 01/24/81 
at 21:16 GMT, (e) 01/25/81 at 09:16 GMT, (f) 01/25/81 at 21:17 
GMT. 

- 200mb wind vectors and vertical motion in pressure coordinate in 
10 mbs (sinking motions are represented by dashed lines and upward 
motion, by sol lid lines) on (a) 01/23/81 at 00:GMT (Vmax = 66 m/s), 
(b) 01/23/81 at 12:GMT (Vmax = 71 m/s), (c) 01/24/81 at 00:GMT 
(Vmax = 77 m/s), (d) 01/24/81 at 12:GMT (Vmax = 72 m/s), (e) 
01/25/81 at 00: GMT (Vmax = 75 m/s), (f) 01/25/81 at 12:GMT 
(Vmax = 68 m/s). 

- Time dependence of heat source. Tmax = 24 hours in the model 
computations. 

- Wind flow at 250mb after 24 hours for a heat source centered at 
10°S and 60°W with symmetrical e-folding width of 1500 km. 

Maximum vector corresponds to 1.98ms” 1 . 

- Same as Fig. 5 but at 48 hours. Vmax = 4.94ms” 1 . 

_ 4 _1 

- Vertical p-velocity in 10 mbs at 500mb after 24 hours. 

- Same as Fig. 7 but after 48 hours. 

- Same as Fig. 7 but for a heat source centered at 15°S and' zonal 
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e-folding width of 700 km and meridional width of 1500 km. 

Fig. 10 - Same as Fig. 5 but at 48 hours for realistic initial flow. 

Vmax = 19.82 ms -1 . 

Fig. 11 - Difference between 250mb wind field at 48 hours and initial flow. 
Vmax - 4.96 ms" 1 . 

- 4 .1 

Fig. 12 - Vertical p-velocity in 10 mbs at 48 hours for a heat source 
located at the Equator and 130°W. 
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